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An  analytical  model  of  penetration  by  non-deforming  pointy  projectiles  into  dry 
sand  is  proposed.  The  model  postulates  that  an  oblique  compaction  wave  in  the 
sand  is  attached  to  the  penetrator’s  tip.  Sand  is  modeled  by  the  P-a  porosity 
model  up  to  full  compaction,  then  is  taken  as  incompressible.  The  sand’s  com¬ 
paction  is  considered  the  most  important  effect;  strength  and  friction  are  ne¬ 
glected.  The  model  handles  wedges,  for  which  a  closed-form  solution  is  derived, 
and  cones,  which  require  numerical  integration.  For  a  given  penetration  veloc¬ 
ity,  the  model  predicts  the  shock  angle  and  the  pressure  on  the  penetrator  sur¬ 
face.  Hydrocode  computations  agree  closely  with  the  model. 

INTRODUCTION 

An  analytical  model  of  penetration  in  dry  sand  by  sharp  conical-nosed  penetrators  is 
proposed.  The  model  predicates  the  existence  in  the  sand  of  an  oblique  compaction  shock 
attached  to  the  nose  of  the  penetrator,  which  is  pointy:  a  wedge  in  plane  flow  or  a  cone  in  ax- 
i symmetric  flow.  The  penetrator’s  deceleration  is  assumed  too  slow  to  affect  the  flow,  so 
from  its  point  of  view  the  sand  flows  steadily  around  it.  Sand,  with  its  intergranular  voids,  is 
modeled  by  the  P-a  compaction  model[l];  for  a  shock  to  form,  the  flow  must  be  supersonic. 

A  compaction  shock  stands  as  an  oblique  wave  attached  to  the  penetrator  apex.  The 
spaces  between  the  sand  grains  are  assumed  to  be  at  least  partially  crushed  out  in  this  wave. 
For  wedge  penetrators,  the  sand  is  deflected  within  the  shock  by  an  angle  equal  to  the  wedge 
angle;  that  is,  the  sand  flows  parallel  to  the  wedge  surface  (and  with  no  further  compaction). 
For  cones,  the  shock  deflects  the  sand  by  an  angle  smaller  than  the  apex  angle,  and  the  sub¬ 
sequent  flow  follows  a  curved  path  that  asymptotically  approaches  the  cone  surface,  with  fur¬ 
ther  compaction. 


MODEL  DEVELOPMENT 


The  compaction  behavior  of  the  sand  is  described  by  the  P-a  model,  Fig.  1,  given  by 

-\2 

P  -P 


a  =  < 


1  +  («o  - 1) 


for  P  <  Ps 
for  P  >  P 


(1) 


where  <2  is  a  porosity  parameter  (with  initial  value  «o  >  1),  a  =  pjp,  in  which  p  is  the  current 
density,  ps  the  fully  compacted  density,  and  Ps  the  pressure  at  complete  compaction.  Parame¬ 
ters  established  for  dry  sand[2]  are  used:  as  pressure  increases,  the  sand  is  gradually  crushed 
from  an  initial  density  of  1 .6  g/cm  to  2.65  g/cm  (ao  =  1 .66)  at  a  pressure  Ps  of  250  MPa. 
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Fully  compacted  material  (quartz) 

P-alpha  model  compaction  curve 

Fig.  1.  The  P-a  model  of  sand  compaction. 


Once  compaction  is  complete,  the  material  is 
assumed  in  this  model  to  behave  incompressibly 
(but  in  the  hydrocode  computations,  fully  com¬ 
pacted  sand  is  modeled  as  quartz). 

Suddenly  compressing  a  porous  material 
creates  a  compaction  wave,  behind  which  some 
or  all  of  the  porosity  has  been  crushed  out, 
depending  on  the  loading.  The  wave  is  actually  a 
shock,  across  which  the  material  incurs  sudden 
changes  in  velocity,  density,  and  pressure,  as 
described  by  the  same  equations  as  govern 
shocks.  The  only  difference  from  the  familiar 
shocks  in  non-porous  materials  is  the 
compaction  behavior  reflected  in  the  porous 
material’s  equation  of  state. 

A  pointy  body  thrust  rapidly  into  a  porous  material  drives  the  material  aside  and  com¬ 
pacts  it  by  a  similar  wave,  which  emanates  at  some  angle  from  the  penetrator  apex.  The  angle 
depends  on  the  apex  angle,  on  whether  the  penetrator  is  a  cone  or  wedge,  and  on  the  penetra¬ 
tion  velocity. 

Relative  to  the  penetrator,  target  material  ap¬ 
proaches  at  the  penetration  velocity  U.  This  veloc¬ 
ity  is  denoted  V\  and  the  velocity  behind  the  shock 
V2.  The  situation  constitutes  supersonic  flow  about 
an  infinitely  long  wedge,  as  shown  in  Fig.  2.  A  par¬ 
ticle  approaching  at  velocity  V\  abruptly  changes  in 
velocity  to  V2  in  a  direction  parallel  to  the  wedge’s 
surface;  thus  it  is  deflected  by  an  angle  equal  to  the 
wedge’s  half- apex  angle  8.  The  shock  makes  an 
angle  cr  with  the  flow  direction. 

The  situation  differs  for  a  cone,  Fig.  3.  The  an¬ 
gle  fiby  which  the  shock  deflects  the  streamline  is 

smaller  than  the  apex  angle  /?,  and  the  flow  behind  the  shock  is  nonunifonn,  as  the  streamline 
curves  to  asymptotically  approach  the  cone  angle  and  the  sand  undergoes  further  compaction. 


Fig.  2.  Oblique  compaction  shock  at  a  wedge. 


Compaction  Wave 

The  development  and  solution  of  the  equations 
follows  standard  texts  on  gas  dynamics  (e.g.,  [3]). 
Changes  in  velocity,  density,  and  pressure  through 
the  shock  can  be  analyzed  with  the  aid  of  Fig.  4, 
which  shows  velocity  vectors  before  and  behind  the 
wave.  Velocity  components  normal  and  transverse 
to  the  wave  front  are  denoted  by  subscripts  n  and  t. 

The  shock  equations  are: 

Continuity  in  //-direction:  pil/tn  =  p2V2n 

Momentum  in  /-direction:  ptVtnVu  =  p2 V2n V2t 
Momentum  in  n -direction:  Pl  +  p{V^  =  P2  +  p2V2n 


Fig.  3.  Oblique  compaction  shock  at  a  cone. 
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In  terms  of  velocity  magnitudes  and  angles,  these  equations  may  be  written  as: 

yOjFj  sin  a  =  p2V2  sin(  cr  -  8) 

pxV2  sin  a  cos  cr  =  p2  V2  sin(cr  -  <7)  cos(cr  -  <7) 

Pl  +  P\V2  sin2  cr  =  P2+  p2V 2  sin2  (cr  -  8) 
to  which  is  added  the  P-a  model,  applied  to  state  2  behind  the  shock, 


■iXi-A/pJ 


for 

for 


Pi<P, 
P2  ^Ps 


(2) 

(3) 

(4) 

(5) 


Equations  (2-5)  relate  cr,  pi  (or  c/2),  Vi,  P2,  and  8. 

For  a  wedge,  the  deflection  angle  8  is  known,  so  a 
solution  may  be  obtained.  For  a  cone,  however,  the 
shock  angle  cr  is  treated  as  given  and  8  as  unknown. 

Squaring  Eq.  (2)  and  dividing  by  Eq.  (3)  yields 
px  tan  a  =  p2  tan(cr  -  8)  (6) 

Solving  Eqs.  (2)  and  (4)  with  Pi  =  0  yields 

pi  =  Pi  sin 2  cr(l  -  px  / p2 )  (7) 

For  partial  compaction  in  the  shock  (P2  <  Ps),  Eqs. 

(5-7)  can  be  combined  into 

tan"’  cr[(a0  -l)(l-5)2  +  l]tan5  +  tan2  a[a0  -  2 B(a0  -1)]+  tan cr[a0  tan 8]  +  a0  =  0  (8) 

where  B  =  p\  V\!P<~.  This  equation,  a  cubic  in  tan  cr,  is  awkward;  further,  it  is  limited  to  par¬ 
tial  compaction  in  the  shock.  It  is  easier  to  assume  a  value  of  cr,  then  solve  for  the  other 
shock  properties;  this  approach  is  sufficient  to  plot  curves.  By  Eqs.  (5)  and  (7), 

— — -=  [l-f?sin2  cr(l-a2/a0)]~ 

« 0  - 1 

which,  given  cr,  may  be  solved  for  the  compaction  parameter  <22, 


the  oblique  compaction  wave. 


(9) 


«1  =  1_ 
an 


1 


B~  sin4  cr 


25  sin 2  cr- 


a„ 


Then 


pi  =  P\V\  sin2  cr(l-«2/a0) 
8  =  cr  -  arctan[(a2  /a0)  tan  cr] 


(10) 

(11) 
(12) 

(13) 

sin(cr  -  8)  a0 

The  above  development  is  based  on  partial  compaction.  As  penetration  velocity  or  shock 
angle  increases,  so  does  the  shock  pressure  P2,  while  cij  decreases,  eventually  reaching  a 
value  of  unity  at  complete  compaction,  for  which  P2  >  Ps,  or,  by  Eq.  (11), 


V2=VX 


sin  cr 


a. 


sin“  cr  > 


a,, 


«0-i  PiK 


(14) 


for  complete  compaction  in  the  shock,  a2=  1 .  Equations  (11-13)  still  hold,  but  Eq.  (6)  re¬ 


duces  to 


which  has  two  roots 


tan  cr  tan  5  -  tancr(a0  - 1)  +  a0  tan5  =  0 


tan  cr  = 


a0  - 1  ±  y/(a0  - 1)2  -  4 a0  tan2  8 


2  tan  8 

If  the  velocity  is  sufficient,  the  curve  of  crvs.  ^intersects  this  curve,  then  runs  along  it. 


(15) 

(16) 
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The  relation  between  shock  and 
deflection  angles  is  shown  in  Fig.  5. 
The  solid  curve  for  complete 
compaction,  Eq.  (15),  has  for  each  8 
two  values  of  cr.  For  8  larger  than 
about  14°,  there  is  no  solution,  as  an 
attached  shock  is  not  possible;  here, 
any  shock  must  be  detached  and 
stand  as  a  bow  shock  ahead  of  the 
penetrator.  Points  on  the  curve’s 
upper  leg  are  not  reachable;  they 
correspond  to  the  so-called  “strong- 
shock”  solution,  but  only  the  lower 
“weak-shock”  value  can  occur. 

The  other  curves  are  partial- 
compaction  curves  for  three  penetra¬ 
tion  velocities.  Of  these,  the  dashed 
curve  corresponds  to  a  critical 
velocity,  as  it  intersects  the  com- 


Complete-compaction  curve 
' ' " '  Partial  compaction,  U  =  4000  ft/s 
Partial  compaction,  U  =  2600  ft/s 
Partial  compaction,  U  =  2000  ft/s 

Fig.  5.  Shock-deflection  solutions  for  three  velocities. 


plete-compaction  curve  at  its  right  end.  This  point  is  located  by  observing  that  Eq.  (16)  has  a 
single  value  here,  so  the  radicand  must  vanish,  which  yields  the  critical  deflection  8C, 

~l  (17) 


tan  &  =  “° 


Substituting  this  into  Eq.  (16)  yields 


tan  cr,  = 


(18) 


For  ao  =  1.655,  these  formulas  give  8C  =  14.3°,  crc  =  52.1°,  which  is  the  rightmost  point  of 
the  complete-compaction  curve.  Thus,  for  wedges  blunter  than  14.3°  half-angle,  an  attached 
shock  cannot  exist,  no  matter  how  high  the  velocity. 

The  penetration  velocity  whose  partial-compaction  o-8 curve  intersects  the  complete- 
compaction  curve  at  this  point  corresponds  to  Pi  =  Ps,  ai=\,  and  cr=  crc,  for  which  Eq.  (11) 

gives  the  velocity  Uc  =  [  5  ao+^  (19) 

VA«  o'1 

At  velocities  above  this,  complete  compaction  occurs  in  the  shock  for  deflection  angles 
greater  than  a  certain  value.  For  the  P-a  model,  this  critical  velocity  is  about  2600  ft/s. 

The  dash-dot  curve  in  Fig.  5  is  the  partial-compaction  curve  for  a  velocity,  2000  ft/s,  be¬ 
low  this  critical  value.  It  intersects  the  complete-compaction  curve  on  its  unreachable  upper 
leg.  At  such  sub-critical  velocities,  complete  compaction  cannot  occur  in  the  shock.  Note  that 
this  curve’s  rightmost  point  lies  left  of  the  rightmost  point  of  the  full-compaction  curve;  for 
deflection  angles  greater  than  this  point’s,  an  attached  shock  is  not  possible,  and  any  shock 
must  be  detached.  The  maximum  deflection  angle  is  given  by  Eq.  (8)  with  d8/da=  0,  or  equi¬ 
valently  by  S( tan d)/S( tana)  =  0,  which  yields  the  condition 

3  tan2  crmax  tan £max [(a0  - 1)(1  -  B)2  + 1]  -  2  tan crmax  [2B(a0  - 1)  -  a0  ]  +  a0  tan £max  =  0 
which  in  conjunction  with  Eq.  (8)  locates  the  point  of  maximum  deflection  for  a  given  veloc¬ 
ity.  Eliminating  tan  from  these  two  equations  yields 
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AC  tan4  crmax  -  tan2  crmax(C  +  3 A)a0  +  «2  =  0 
where  A  =  (or0  - 1)(1  -  B) 2  +1  and  C  =  2B(a0  - 1)  -  a0,  for  which  the  solution  is 


tan"  cr  „  =  a. 


C  +  3A  +  tJ(C  +  3A)2  +  4  AC 


2  AC 


The  maximum  deflection  angle  for  an  attached  shock  is  then  given  by 

2[2 B(a0  -l)-a0]tano-max 


tan  />'  = 


<j  +  an 

max  0 


(20) 


(21) 


4500  “ 


4000  - 


«  3500  - 


i - 1 - n - 

i 

Attached  shock, 
complete  s  180 
compaction 


3[(a0  -l)(l-5)4  +  l]tan' 

For  a  penetration  velocity  of  2000  ft/s,  this  angle  is  about  11°. 

The  dotted  line  in  Fig.  5  is  the  partial-compaction  curve  for  a  velocity,  4000  ft/s,  above 
that  given  by  Eq.  (19).  It  intersects  the  complete-compaction  curve  along  its  lower  leg,  at  a 
deflection  angle  of  about  11°.  For  deflection  angles  smaller  than  this,  only  partial  compaction 
occurs  in  the  shock;  for  larger  deflections,  complete  compaction  occurs  there. 

Each  partial-compaction  curve  in  Fig.  5  meets  the  vertical  axis  at  a  positive  value,  the 
minimum  possible  shock  angle,  which  corresponds  to  an  infinitesimal  deflection  or  a  shock 
of  vanishing  strength  traveling  at  the  sound  speed  c0  =  y[dP~! dp  .  Differentiating  Eq.  (1)  and 
evaluating  at  a  =  «o  yields  Co  =  aoPs  /[2p0  («o  -  1)].  For  a0  =  1 .66,  Ps  =  250  MPa,  p0=  1.6 

n  _ i 

g/cnr,  the  sound  speed  is  443  m/s.  The  minimum  shock  angle  is  the  Mach  angle,  shf  ( co/U ). 

This  model  requires  an  attached  shock,  whose  existence  depends  on  the  velocity  and  de¬ 
flection  angle.  At  velocities  below  the  sound  speed  co,  no  shock  is  possible.  At  velocities 
above  the  critical 
value  Uc  of  Eq.  (19),  5000 

the  shock  is  attached 
for  deflection  angles 
up  to  8C  of  Eq.  (17). 

Between  co  and  Uc, 
the  shock  is  attached 
for  deflection  angles 
up  to  Amax  given  by 
Eq.  (21). 

Curves  for  these 
equations  divide  the 
plane  of  velocity  vs. 
deflection  angle  in 
Fig.  6  and  define  the 
limits  of  applicability 
of  the  attached-shock 
model.  The  model 
holds  at  all  points 
above  the  dashed 
curve,  which  corre¬ 
sponds  to  Eqs.  (20- 
21);  below,  the  shock 
is  detached  or  non¬ 
existent  (in  the  sub- 


Eqs.  (14)  &  (16) 


3000 


2500 


2000  - 


1500 


1000 


Attached  shock, 
partial  compaction 


Eqs.  (20)  &  (21) 


UCW808.4 


Detached  shock  Cn-)2808.4 


_L 


_L 


Subsonic,  no  shock 

J _ I _ L 


_L 

12 


2  4  6  8  10 

Wedge  angle  (degrees) 

Fig.  6.  Regimes  of  shock  behavior  over  the  plane  of 
penetration  velocity  vs.  shock  deflection  angle. 
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sonic  regime  at  the  bottom).  Above  the  solid  curve  of  Eqs.  (14)  and  (16),  the  shock  com¬ 
pletely  compacts  the  sand.  Between  these  curves,  the  shock  only  partially  compacts  the  sand. 

Conical  Flow  Behind  the  Shock 


A  conical  penetrator  differs  from  a  wedge  in  that  the  flow  is  axisymmetric,  and  behind 
the  shock,  as  in  Fig.  3,  each  streamline  curves  as  it  approaches  the  cone.  This  flow  is  deter¬ 
mined  from  the  solution  of  Taylor  and  Maccoll[4]  for  supersonic  flow  about  a  cone.  The 
same  dynamical  equations  are  used  here,  but  the  P-a  model  replaces  the  ideal  gas  law. 

The  solution  is  based  on  the  observation  that,  for  flow  about  a  cone  of  infinite  extent,  the 
problem  involves  no  length  dimension;  the  solution  therefore  depends  only  on  angles  cen¬ 
tered  on  the  cone’s  apex.  Using  a  spherical  coordinate  system  centered  here  reduces  the 
equations  for  the  flow  behind  the  shock  to  the  ordinary  differential  equations, 


2Vr  +  Ve  cot  9  + 


dVa 


V0  da 


dd  a  d6 


=  0;  V6=dVr/d0; 


dP_ 

d6 


Ps_ 

a 


V. 


dVr 

dd 


+  va 


dV g_ 

d6 


(22) 


where  Vr  and  Vq  are  velocity  components  in  the  radial  and  angular  directions  (measured  from 
the  apex).  The  boundary  conditions  are  the  states  behind  the  shock  given  by  the  oblique- 
shock  relations,  and  Vq  =  0  at  the  cone’s  surface. 

The  solution  scheme  is  inverse.  For  a  given  U,  a  value  of  shock  angle  o  is  assumed.  The 
conditions  behind  the  shock  are  given  by  Eqs.  (10-13).  Velocity  components  there  are  given 
by  V<\e=a  =  ^2  sin(<r  -  S) ;  Vr\e=o.  =  V2  cos(cr  -  S) .  Eqs.  (22)  are  then  numerically  integrated 

from  the  shock  to  the  cone  for  decreasing  6,  until  Vq  =  0,  which  satisfies  the  boundary  condi¬ 
tion  at  the  surface,  where  the  value  of  0  is  the  apex  angle  fd  that  yields  this  solution.  At  low 
velocities,  the  sand  is  only  partially  compacted,  but  at  increasing  velocities,  complete  com¬ 
paction  occurs  first  near  the  cone  surface,  then  encompasses  an  increasing  angle,  until  com¬ 
plete  compaction  occurs  in  the  shock. 

Calculated  results  at  3000  ft/s  are  shown  in  Fig.  7,  which  plots  the  relationships  among 
the  various  angles. 

The  dark-blue  curve 
relates  shock  angle 
to  deflection  angle 
at  full  compaction. 

The  green  curve  in¬ 
tersecting  the  verti¬ 
cal  axis  at  the  Mach 
angle  of  29.1°  re¬ 
lates  the  shock  and 
deflection  angles;  it 
starts  to  overwrite 
the  blue  curve  at  a 
shock  angle  of 
43.4°,  where  full 
compaction  in  the 
shock  begins.  Its 
right  end  corre¬ 
sponds  to  the 


eo- 


Shock 

angle, 

a 

(deg) 

40- 


20- 


j  c< 

maximum  d 

implete  comj 

e  flection 

laction  in  she 

complet 

e 

\ 

Mach 

angle 

Full 

compaction 

compac 

boundai 

c 

§  * 

5  10  15  20  25 

Deflection  angle,  8,  or  cone  apex  angle,  ft  (deg) 

Fig.  7.  Computed  results  for  a  penetration  velocity  of  3,000  ft/s. 
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maximum  deflection  for  an  attached  shock,  14.3°,  at  a  shock  angle  of  52.2°.  The  red  curve 
relates  shock  angle  to  cone  apex  angle;  because  the  maximum  shock  angle  is  52.2°  (upper 
dashed  line),  the  uppermost  point  reachable  on  this  curve  is  at  an  apex  angle  of  27.2°.  The 
light-blue  curve  locates  the  boundary  between  partially  and  fully  compacted  sand;  for  apex 
angles  less  than  about  22°,  no  sand  is  fully  compacted,  but  as  //  increases,  this  boundary  ex¬ 
pands  outward  from  the  cone  surface  until  full  compaction  occurs  in  the  shock,  where  it  in¬ 
tersects  the  shock-angle  curve  (lower  dashed  line)  at 
about  a  24°  apex. 

COMPARISON  WITH  HYDROCODE 

To  compare  with  the  model,  the  CTH  hydro- 
code[5]  was  run  for  cones  penetrating  dry  sand 
modeled  by  the  P-a  model.  Figure  8  shows  the  re¬ 
sult  for  a  20°  cone  at  3,000  ft/s.  The  plot,  of  color 
contours  of  density,  reveals  an  oblique  wave  (in 
green)  attached  to  the  apex  of  the  (red)  cone.  The 
model  predicts  a  shock  angle  of  about  36. 1°,  in  good 
agreement  with  the  angle  in  this  plot  (35.5°). 

Computations  were  run  for  apex  angles  of  10° 
to  30°  and  velocities  of  3,000  and  4,000  ft/s.  Values 
of  shock  angle  cr  measured  from  the  CTH  plots  are 
compared  with  model  predictions  in  Table  1;  agree¬ 
ment  is  within  a  fraction  of  a  degree  in  all  but  one 
case.  For  the  30°  cone,  the  model  predicts  detached 
shocks  at  both  velocities,  and  the  CTH  plots,  as  in 
Fig.  9,  clearly  show  curved  waves.  The  CTH  values 
are  shown  on  the  model  plots  in  Fig.  10. 


Fig.  8.  CTH  computation  of  a  20° 
cone  penetrating  sand  at  3,000  ft/s. 


Table  1.  Comparison  of  analytical  model  and 
CTH  computations  for  conical  penetrators. 
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Fig.  9.  CTH  computation  of  a  30°-half-angle 
cone  penetrating  dry  sand  at  3,000  ft/s. 


CONCLUSIONS 


An  analytical  model  for  pointy  wedges  and  cones  penetrating  dry  sand  at  supersonic  ve¬ 
locities  has  been  developed.  The  model  predicts  the  pressure  acting  on  the  penetrator. 
Agreement  with  hydrocode  computations  is  excellent. 
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Fig.  10.  Model  solutions  (curves)  and  CTF1  predictions  (symbols)  of  shock  angle  for  penetration 

at  3,000  ft/s  (top)  and  4,000  ft/s  (bottom). 
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